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ABSTRACT 

The relationship between the metric and nonrelativistic matter distribution depends 
on the theory of gravity and additional fields, hence providing a possible way of distin- 
guishing competing theories. With the assumption that the geometry and kinematics of 
the homogeneous Universe have been measured to sufficient accuracy, we present a pro- 
cedure for understanding and testing the relationship between the cosmological matter 
distribution and metric perturbations (along with their respective evolution) using the 
ratio of the physical size of the perturbation to the size of the horizon as our small ex- 
pansion parameter. We expand around Newtonian gravity on linear, subhorizon scales 
with coefficient functions in front of the expansion parameter. Our framework relies 
on an ansatz which ensures that (i) the Poisson equation is recovered on small scales 
(ii) the metric variables (and any additional fields) are generated and supported by 
the nonrelativistic matter overdensity. The scales for which our framework is intended 
are small enough so that cosmic variance does not significantly limit the accuracy of 
the measurements and large enough to avoid complications due to nonlinear effects 
and baryon cooling. From a theoretical perspective, the coefficient functions provide a 
general framework for contrasting the consequences of ACDM (cosmological constant 
-|- cold dark matter) and its alternatives. We calculate the coefficient functions for 
general relativity (GR) with a cosmological constant and dark matter, GR with dark 
matter and quintessence, scalar-tensor theories (STT), f{R) gravity and braneworld 
(DGP) models. We identify a possibly unique signature of braneworld models. For 
observers, constraining the coefficient functions provides a streamlined approach for 
testing gravity in a scale dependent manner. We briefly discuss the observations best 
suited for an application of our framework. 
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1 INTRODUCTION 

A successful model of the Universe must include a back- 
ground geometry, an inventory of its contents, a kinematical 
description of its expansion and a dynamical explanation 
of liow its constituents interact, drive the expansion and 



develop structure. Recent observations (for example 


Riess 


et al. |1998 


Perlmutter et al. |1999 Preedman et al.| 


2001 


Allen et al. 2002"Tegmark et al."2004"Riess et al.'2004 
et al. 2005; Eiscnstein et al.;2005; Astier et al. 2006; Sj 


'Cole 
jergel 


et al.| 2007 


Percival et al.| 2007, Allen et al. 


|2008| IMantz 


et al.|2008| 


Komatsu et al.|2008jand references therein) have 



led to a "Flat ACDM cosmology" (henceforth FACDM), 
dominated by dark energy (cosmological constant A) and 
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matter (predominately dark and initially cold) and the ob- 
served expansion rate and growth of structure agree with 
the predictions of this model at the ten percent level. Fu- 
ture observations should be capable of testing this model 
at the one percent level. If they verify its predictions, they 
will affirm a remarkable, simple description of the Universe, 
implicit in the earliest relativistic investigations of Einstein, 
Friedmann and Lemaitre, analogous to the affirmation of 
general relativity (GR) that took place twenty years ago 
(for example [Win||2006| ). If FACDM passes this test, then 
the challenge will be to account for this outcome in terms 
of physical processes operating at earlier epochs; if it fails, 
then we shall either have learned something important about 
gravitational physics or described a new, dominant compo- 
nent of the Universe. Many alternatives, with and without 
GR, to FACDM have been proposed. At this stage, none 
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of them stands out. There is therefore a need to provide a 
framework for describing future observations and theoretical 
investigations in general terms which will facilitate a distinc- 
tion between FACDM and its alternatives. The provision of 
one such framework is the goal of this paper. 

Further observational progress is anticipated over the 



coming decade. The analysis of Planck observations | Planck- 
Collaboration||2006| ) of the microwave background, coupled 
with local measurements of the contemporary Hubble pa- 
rameter, Ho, should result in an extremely accurate descrip- 
tion of the physical conditions and the statistical properties 
of the density fluctuation spectrum at the epoch of recombi- 
nation when the Universe had a scale factor a= {l + z)~^ ~ 
10~^ relative to today. Combining the calculated physical 
sizes of the acoustic peaks in the background radiation spec- 
trum with the Hubble constant and the Copernican Princi- 
ple leads to a measurement of spatial curvature, which is al- 
ready known to contribute to the kinematics at a level of less 
than a few percent ( Komatsu et al.|2008 l. We shall adopt a 
value of zero for illustration purposes. Essentially kinematic 
measurements, for example, those involving Type la super- 
nova explosions, baryonic acoustic oscillations (BAO) and 
baryonic gas fractions in clusters should provide a record 
of the comoving distance, d{a) = J cdt/a, from which the 
evolution of the Hubble parameter H{a) = dlna/dt and 
the acceleration parameter q{a) — d In(Ha) /d In a can be in- 
ferrecj^ For the rest of the paper we shall assume that these 
evolutions have been measured to a sufficient accuracy. Note 
that we are using a instead of the cosmic time t as the time 
coordinate as this relates directly to the observable photon 
frequency shift. For recent constraints on the expansion his- 
tory, see for example (Rapetti et al 
therein. 

Given an understanding of the geometry and kinemat- 
ics, the task is then to see if the dynamical evolution of 
the Universe is consistent with GR or mandates an alterna- 
tive theory. Now, GR provides a relationship between the 
spacetime geometry on a cosmological scale measured by 
the Einstein tensor G[(?^,y] and the total Energy-Momentum 
Tensor (EMT) of its contents T, G[ff^^] = SttGT. The dis- 
covery that G[g^„] / 87rGT["obs"] where T["obs"] includes 
known forms of matter such as electromagnetic radiation, 
baryons etc. has led to the addition of dark matter and dark 
energy contributions to the EMT. Dark matter candidates 
include Weakly Interacting Massive Particles and axions 
which would presumably behave gravitationally like bary- 
onic matter. However other possibilities exist which might 
behave differently such as massive neutrinos (as a subdomi- 
nant component). Dark energy is most simply characterized 
as a temporally and spatially constant vacuum energy field 
with zero enthalpy (see ICarroU 2001 for a review). However, 
it could also have quite difi'erent dynamical properties and 



20071 and references 



might include contributions from additional scalar (Ratra 
fc Peebles|1988[ ), vector ( |Armendariz-Picon|2004[ ) or tensor 
fields with possible interactions between each other (Farrar 
fc Peebles||2004 \ and with known forms of matter. Histori- 
cally, the first representation of dark energy was Einstein's 
cosmological constant, which was seen as an augmentation 



^ Our acceleration parameter differs from the conventional decel- 
eration parameter by a minus sign. 



to G, not T (see for example Carroll et al.|20d4| . This origi- 
nal proposal has also been generalized in many ways so that 
G[g^v]+F[g^^, (p] = 87rGT["obs"], where F[g^^, tp] depends 
on the metric and more generally some additional gravita- 
tional fields, if. For example <p could be the additional gravi- 
tational scalar field in Scalar- Tensor Theories (STT) (see for 
example Santiago et al.|[l998 Perrotta et al.pOOO l. Nature 
could of course be unkind and we might have 



G[gf,^] + F[g^,, <p] = 87rGT["obs"] + SttGT ["dark" 



(1) 



Considerable effort has been made in constructing models 
that fall into the above mentioned categories and more re- 
cently in finding ways to distinguish between them (for ex- 
ample see|Lue et al.|2004||Ish ak et al.'2006', [Bludman|2007l 
IZhang et al.||2007l IHuterer & Linder 2007). 



Now, modifying the physics beyond GR with cold dark 
matter and A can have three quite separate manifestations. 
Firstly it can lead to a change in expansion of the universe, 
secondly, it can infiuence the growth of structure and the 
metric and thirdly, it can confront local tests of the theory 
of gravity. The approach that we follow is to assume that the 
theory is constrained by the first and third manifestations 
and that it is the growth of structure that is providing the 
test. This oversimplifies the data analysis but does lead to 
a transparent and simple approach. One important conse- 
quence of adopting local gravitational tests is that photons 
and baryons, at least, will follow geodesies and that the un- 
perturbed photons will be subject to cosmological redshift- 



ing of their frequencies. 



This simplifies the inter- 



pretation of observational data. 

Our procedure is to adopt a general form for the metric 
of a linearly perturbed homogeneous and isotropic universe 
which introduces two potentials "I>(x, a) and 5'(x, a) (scalar 
metric perturbations in the Newtonian gauge), where x de- 
notes the three spatial coordinates. We also introduce an 
associated fractional density perturbation 5m (x, a) in non- 
relativistic matter and relate it to the potentials. We assume 
that there is a dominant nonbaryonic contribution to the 
clustering of nonrelativistic matter. In practice, it is easier 
to work with Fourier modes and this allows us to focus at- 
tention on the range of length scales that are most relevant 
observationally: sufficiently smaller than the horizon so that 
our expansion is valid and we can observe enough indepen- 
dent volumes within our current horizon allowing for a high 
precision measurement despite "cosmic variance", and yet 
large enough that nonlinear effects and baryonic cooling are 
not a factor. Within this range of length scales, we adopt 
the following ansatz regarding the relationship between lin- 
earized metric and density perturbations, written as an ex- 
pansion in powers of (aH/k), where k is the magnitude of 
the comoving wavevector k 



$(k,a) 



<5m(k, a)B 



aH 
aH 



(2) 



(5™(k, a) = (5mi(k)A 



aH 
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where 



B 


aH 
















aH 




r 














'aH 




A 













/3o(a) + /3i(a) ( ^"j + p^ia) +, 



such as effective gravitational constant on small scales ( Tsu- 
jikawa||2007[), growth of perturbations on s mall scales (|Lin- 



70(a) +71 (a) + 72(a) 



+ . 



5o{a) + 5,{a)['^]+5,{a)[^-f\ +. 



aH 
k 



We have set the speed of light c = 1. The background mass 
density pm. oc a'"^. 5mi{k) is determined from initial condi- 
tions which can in principle be taken close to the surface of 
last scattering, ~ 10""^, as long as the modes are suffi- 
ciently sub-horizon. Note that often, we are only interested 
in the scale dependence of the growth of the perturbations 
in the matter distribution and the metric at linear, sub- 
horizon scales. Such measurements require taking ratios of 
the matter distribution or the metric at different redshifts, 
whereby, the initial conditions cancel out in the final expres- 
sions. The coefficient functions {/3n,7n,<5n} with n = 0, 1,2 
are arbitrary functions of the scale factor. The leading terms 
in the expansion agree with Poisson's equation on small 
scales, while subsequent terms allow for a scale-dependent 
departure as we move towards larger scales. This approach 
introduces a length-scale dependence to the perturbations 
through an expansion in powers of (GM/dp)^^^ ~ dp/dn, 
where M ~ pmdp is the total mass enclosed within the phys- 
ical size dp and du = 1 /H is the Hubble horizon. In Fourier 
space, with dp ~ a/k, we get (GM/dp)^^^ ~ aH/k. 

The theories that we discuss below introduce differ- 
ent corrections (different coefficient functions {/3n,7n,5n}) 
and these differences are measurabl^E] • From an observer's 
perspective, constraining the coefficient functions with mea- 
surents of "1>, and 5m provides a streamlined approach to 
characterizing gravity on cosmological scales in a scale de- 
pendent manner. On the other hand, from a theorist's per- 
spective, substituting the ansatz into the field equations for 
a given theory allows for a (mostly straightforward) calcu- 
lation of the coefficient functions. The coefficient functions 
provide a means of comparing the consequences of different 
theories. We shall discuss our assumptions, limitations and 
our ansatz in detail in the next section. 

This is certainly not the first time that an attempt 
at constructing and applying such a framework has been 
made. The Parametrized Post Newtonian formalism (see 
WillpOOSl and references therein) has been a powerful frame- 
work for understanding and constraining gravity on solar 
system (and other isolated system) scales. Our aim is to con- 
struct a similar framework for cosmological scales. Recently 
a few attempts have been made in this direction. However 
most of these are either concerned with the expansion his- 
tory alone, deal with specific aspects of departures from GR 



^ As we shall see at the end of Section 2 and in Section 3, the 
coefficients /32 , 72 1 <52 7^ even in GR with non-relativistic mat- 
ter and represent "post-Newtonian" corrections. Also note that 
Sm characterizes the fractional matter density perturbation in 
the Newtonian gauge, which is related to the often used comov- 
ing density perturbation Am through Am = &m + 3ada{5m — 
3'^){aH/k)^ . A combination of the 00 and Oi Einstein equations 
yields (k/aHf^ = -{iTTGpm/ H^)Am- 



der & Cahn 20071, the gravitational slip (Caldwell et al 
2007 1 , or deal with superhorizon scales ( Bertschinger|200"6 1 



The authors in ( [Amendola e t aC 2008) take into account 
growth of structure, anisotropic stress and the modification 
to the Poisson equation and parametrize departures from 
Einstein's gravity with a growth index and two functions 
of the scale factor which are relevant for weak lensing sur- 
veys. However, they do not consider scale dependent de- 
partures. Another popular phenomenological approach for 
characterizing the effects of the unknown physics (additional 
fields, their interactions, or modified gravitational laws) is 
to define an effective fluid energy momentum tensor for ev- 
erything other than the standard model matter, effectively 
move F in equation ([T]) to the right hand side and define 
Te// = -(87rG)~^F[v3,g^^] + T["dark"]. This effective en- 
ergy momentum tensor is then parametrized in terms of the 



equation of state, sound speed, anisotropic stress, etc. (Hu & 
|Eisensteiii||1999 Bashinsky||2007 1 . This approach, however. 



seems to put an unnecessary restriction of a fluid interpreta- 
tion which might be misleading, especially when the effective 
dark energy is due to modifled gravity or extra dimensions. 
We are unaware of a systematic approach undertaken where 
the framework includes a scale dependent departure in the 
relationship between the matter distribution and the met- 
ric perturbations along with their respective evolution on 
cosmological scales up to post-Newtonian order. |^ 

The rest of the paper is organized as follows. Section 
2 discusses our assumptions and the particular form of the 
ansatz in detail. In Section 3 we apply our framework to GR, 
STT, q uintessence, f{R) models f^axroU 2001) and DGP 
gravity (Dvali et al. 20001. In particular, we calculate the 



coefficient functions in these theories and comment on our 
ansatz in the context of these theories. Section 4 is devoted 
to how our framework might be employed by observers. We 
briefly discuss the observations that could be used to con- 
strain the different coefficient functions. Section 5 presents 
a short summary and future directions for extending the 
framework. 



2 OUR ANSATZ AND ASSOCIATED 
ASSUMPTIONS 

With an eye towards observations in the next decade, we 
assume that the geometry (spatial curvature) and kinemat- 
ics (expansion history) of the Universe have been measured 

We note that during the final stages of preparation of this this 
paper we became aware of a scale dependent framework for mod- 



ified gravity that includes super and sub-horizon scales | |Hu fc| 
|Sawicki|2007[ l. After submission of this manuscript, the following 
were posted on arXiv.org which are relevant to this work. | |Jain| 
I&: Zhang|2007[ l provide an analysis of the observational tests for 
modified gravity; l |Hui fc Parfrey||2007| use evolution of galaxy 
bias to constrain scale dependent departures from GR; whereas 
( Bcrtschin ger fc Zukin|2008[l build on jBertschi ngcr 2006 ~| to in- 
cludc sub-horizon scales; | |Hu|2008l l extend l |Hu fc Sawicki|2007[ l 
to include multiple fiuids and curvature relevant for cosmic mi- 
crowave background calculations and constraints; whereas | |Cald-| 
[well et al.||2007l l discuss the effects of gravitational slip on the 
CMB, growth of structure, and lensing observations. 
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to a percent level accuracy. What remains to be understood 
and measured accurately (at the few percent level) is the 
relationship between the metric fluctuations and the non- 
relativistic matter distribution along with their respective 
evolution on linear, subhorizon scales. This relationship will 
depend on the theory of gravity or the presence of yet un- 
known components, thus providing a test for distinguish- 
ing different theories. To explore this relationship in an (al- 
most) model independent way, we provide an ansatz, equa- 
tion relating the scalar metric perturbations (in New- 
tonian gauge) and the nonrelativistic matter overdensity in 
Fourier space. In this section we discuss the particular form 
of the ansatz and the underlying assumptions in detail. We 
introduce our notation and conventions followed by some 
physical arguments regarding our choice of the particular 
form of the ansatz. We end with a discussion of the range of 
scales for which our ansatz is expected to be useful. 

We focus on a perturbed FRW universe (spatially flat) 
with scalar metric fluctuations in the Newtonian gauge 
| Bardeen|1980 1. In this gauge the metric takes the following 
form (c = 1) 

ds^ = -[1 + 2^{yL,t)]dt^ + {t)[l - 2*(x, f)]dx ■ dx 

Here the metric perturbations |<l>(x, t)|, [^^(x, t)| <C 1. We 
choose to work in the Newtonian gauge because $(x, t) is 
the generalization of the Newtonian gravitational poten- 
tial and the potentials $(x, t) and 'I'(x, t) are gauge in- 
variant Bardeen variables when we specialize to the New- 
tonian gauge. The energy density perturbation (5m(x, t) is 
also gauge invariant, corresponding to the energy density 
perturbation on the zero shear spatial hypersurface which 
is closest to Newtonian time slicing [see equation (3.14) in 
(Bardeen 1980 1]. In what follows, we use the scale factor 
a as the independent variable instead of cosmic time t with 
a(today) = 1. With this change of variables, the metric takes 
the form 

ds^ = "[1 + 2$(x, a)\{aHy'^da + a[l ~ 2*(x, a)]dx ■ dx 

We shall work primarily in Fourier space and use the con- 
vention f{x,a) = (27r)~^/ d^k/(k, a)e''' ''. To avoid unnec- 
essary clutter we write the Fourier transform of the metric 
perturbations "^(k, a)e'^'^ as $. The same is true for 4' and 
S,n. The background quantities depend on a. We shall often 
suppress this dependence; for example by H we mean H{a). 

We have assumed spatial flatness as expected on the ba- 
sis of the simplest interpretation of inflation. If the Universe 
has measureable spatial curvature or large scale deviations 
from the Robertson- Walker assumptions of homogeneity and 
isotropy, then the following development must be general- 
ized at the expense of introducing parameters that need flt- 
ting. A purely geometrical demonstration of spatial flatness 
would obviate some of this concern. Such a demonstration is 
possible, in principle, using two screen gravitational lenses 
(Blandford 2008, in preparation), though it is not known 
how practical it will be to implement this demonstration. If 
we choose to include curvature as an additional parameter, 
then location of the flrst acoustic peak in the CMB (and 
BAO scale) would likely provide the best constraints. 

Our ansatz provides a relationship between $, 'I' and 
5m on linear (in <1>, 5' and (5m), subhorizon scales. We now 
turn to the discussion of some important features of this 
ansatz. On scales that are much smaller than the size of the 



horizon, aH/k <^ 1, the leading term has the form of a lin- 
earised Newtonian gravitational fleld equation. For the pur- 
pose of this paper the Newtonian form of the field equation 
refers to the the following relation between the time-time 
metric perturbation "l>(x,a) and the nonrelativistic matter 
density contrast 5m(x, a), V^"I'(x, a) oc 5„i(x, a), which in 
Fourier space becomes <& oc {aH/k)^5m- Now, in the New- 
tonian gauge tl>(x, a) plays the role of the Newtonian po- 
tential once the background has been subtracted out. The 
proportionality allows for a possible temporal variation in 
the effective Newton's constant which could depend on the 
cosmological background evolution. 

From GR we know that this Newtonian relation starts 
breaking down as the size of the perturbation becomes com- 
parable to the size of the horizon. In general, different the- 
ories of gravity will introduce different scale dependent de- 
partures from this equation, changing the metric-matter re- 
lationship. Our claim is that for a large class of theories, 
our ansatz, equation ([2|, captures the scale dependence of 
the relationship between the nonrelativistic matter distribu- 
tion and cosmological metric perturbations. In particular, 
our ansatz faithfully reproduces the scale dependence of the 
metric-matter relationship in the fiducial case of GR with 
cold dark matter and a cosmological constant. In the pres- 
ence of additional fields one might expect this relationship 
to break down; however, this is usually not the case. Suppose 
that an additional field enters the equations, for example as 
a source (quintessence), as a time varying gravitational con- 
stant (Brans Dicke theory) or indirectly encapsulating the 
effect of higher dimensions, etc. Perturbations Sip in such a 
scalar field (p (consider quintessence or scalar-tensor theo- 
ries) will be involved in the relationship between Sm and "1?. 



How ever, from the field equation for Sip, equations ( 18 1 and 
(19 1, we can see that 5(p oc ^{aH/k)^ for quintessence and 



5iy5 oc $ for scalar-tensor theories when aH/k <^ 1. Thus, 
even if additional scalar fields are present, our ansatz should 
be a good approximation for the relationship between the 
matter distribution and the metric at the scales of interest. 
Note that we have assumed 'i' — 0[^] for this argument. 

Another feature of our ansatz is that "1? and ^' are di- 
rectly proportional to 5m- This might seem unusual, since 
it implies that in the absence of nonrelativistic matter per- 
turbations, there would be no metric perturbations. This is 
certainly not true in principle if an additional scalar field is 
present. However observationally, we know that nonrelativis- 
tic perturbations are present and they dominate over per- 
turbations in other fields. The following argument provides 
a more detailed justification. Since on the smallest scales, to 
lowest order in {aH/k), the potential $ oc 5m{aH / k)^ , we 
have Stp oc 5m{aH / k)'*' and Sip oc Sm{aH/k)^ in quintessence 
models and STT respectively. This means that the poten- 
tials and pertubations in other scalar fields are supported 
by the nonrelativistic matter perturbations. We do not ex- 
pect to see the effects of the initial power spectrum of these 
additional fields up to the order of the terms considered in 
our ansatz, with the initial power spectrum of the additional 
field possibly playing a role in higher order terms. This is one 
of the reasons for not extending the power series in aH/k 
beyond the order considered in the ansatz. 

Our ansatz does not capture the matter-metric relation- 
ship for all available models in the literature. Consider for 
example fc-essence ( [Armendariz-Picon et al.|2000[ ), where the 
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effective "sound speed" (cs) can be small. This leads to a sig- 
nificant clustering of dark energy on small scales which can 
be comparable to nonrelativistic matter perturbations. In 
these scenarios, our ansatz does not provide a good approxi- 
mation to the full theory. The coefficients /32 , 72(oc cf^) ^ 1 
signaling a breakdown in our assumptions. More generally, 
if a model introduces an additional physical scale within 
the range of scales of interest, then care needs to be taken 
in using our ansatz. In the fc-essence example, this addi- 
tional scale is the Jean's length for the scalar field fluctu- 
ations, whereas in the case of f{R) models this could the 
"Compton wavelength" (~ /jj^. ) of the effective gravita- 
tional scalar degree of freedom (see for example Pogosian & 
|Silvestri^2008J . In such cases the ansatz might still be ap- 
plicable in a more limited range of parameters and length 
scales (see Section 3.4). 

We note that some of the above arguments are made 
under the assumption that the additional gravitational or 
nongravitational contribution to the field equations is due 
to a scalar field (quintessence or scalar tensor theories). As 
argued above, this leads to only even powers of aH/k in the 
expansion. Furthermore, f{R) modiflcation of the Einstein- 
Hilbert action also lead to even powers of aH/k. An intrigu- 
ing case where one can get an odd power of aH/k is in DGP 
braneworld models. In these extra-dimensional theories, the 
junction conditions on our 4 dimensional brane gives rise to 
a scale dependence involving terms linear in aH/k. We come 
back to this in Section 3.5. 

Regarding we assume that the relationship between 
^ and 5m has the same {aH/k) dependence as $ and 5rn 
since from GR we expect <1> = ^' when no anisotropic stress is 
present. The form of 5m (k, a) in the ansatz can be motivated 
from the conservation equation for nonrelativistic matter at 
first order in and Srn- 



999 315.2 99 





Figure 1. The ratio of the physical size of the perturbation to 
the size of the horizon is used as an expansion parameter in our 
anzatz. We plot this ratio, {aH/k), as a function of a from last 
scattering to the present for the concordance model (yellow re- 
gion) . The upper and lower bounds of the yellow region are deter- 
mined by considering scales that are small enough so that cosmic 
variance does not dominate the errors and at the same time large 
enough so that nonlinear evolution and baryon cooling are not 
a significant factor. Most of the observations in the next decade 
will yield information in the range 10~^ < a < 1. If we are in- 
terested in observations that only care about a smaller range of 
the scale factor, then the allowed range of Ha/k increases. We 
also plot lines of constant multipole I ^ kd(a), which provides a 
rough estimate of the relationship between k and angular scales 
at different redshifts. 



a daSm + {2 + q)adaS„ 



k 
aH 



(3) 



+ 3[a^dl'i> + i2 + q)ada'i'] 



As discussed above at lowest order in (aH/k)^, the metric 
perturbations "I>, 4' oc Sm{aH/k)^ , thus the largest term on 
the RHS of equation ([sjl is proportional to 5m (k, a). At this 
order we get a homogeneous equation for 5m which has a so- 
lution of the form (5m(k, a) — 5mi(k-)5o{a). This is the usual 
approximation used when investigating the growth function 
on small scales. Perturbatively including the next order term 
on the RHS, we can see that our ansatz captures the general 
form of the solution to that order. Again, we use this argu- 
ment as motivation for the form of the ansatz, being aware 
of the fact that nonrelativistic dark matter is not covariantly 
conserved in some models. In 5m, we include both baryonic 
and nonbaryonic dark matter, with an understanding that 
baryonic matter contibutes a small fraction to the total. We 
assume that baryons are covariantly conserved and follow 
timelike geodesies, serving as test particles whose motion 
can be used to probe the metric. 

We now turn to a discussion of the range of scales where 
we expect our procedure to be applicable. Our ansatz uses 
the ratio of the physical size of the perturbation dp{a) to 
the size of the Hubble horizon dH{a) = 1/ H(a) as our small 
(post-Newtonian) expansion parameter. In Fourier space 
dp (a) ~ a/k and we need dp{a) / dni^a) ~ aH/k <^ 1 for 



the expansion in aH/k to be meaningful. From Figure 1, 
we can see that for a given k, aH/k is a decreasing func- 
tion of the scale factor (till dark energy domination) . So if 
aH/k ^ 1 at early times, it will remain so till today. 

We flrst give a rough upper bound on Ho/k. In addition 
to aH/k <^ 1, from an observational standpoint, the largest 
scales of interest are the ones where cosmic variance does not 
significantly limit the precision of our measurements (angu- 
lar multipoles I ^ 1). To convert this constraint on / to a 
constraint on Ho/k, we need a relationship between k and 
I. For large k, a perturbation with a given k corresponds 
roughly to a multipole I ~ kd{a) — {aH/k)~^aHd{a), where 
d{a) is the co-moving distance. Note that this relationship 
is scale factor dependent. Let us take / ~ 30 as the largest 
angular scale where cosmic variance does not signiflcantly 
limit measurement precision. For 0.1 < a < 0.5 we get 
aH/k ~ aHd{a)/l < 0.06 since aHd{a) < 3.7 in this 
range. At a ~ 0.5, the corresponding comoving wavevec- 
tor of the perturbation is fc ~ 10~^^Mpc~^ or equiva- 
lently Ho/k ~ 3 x 10~^. On the other hand, this same I 



would correspond to aH/k > 1 for a 



10- 



So if we 



are also interested in the CMB, then aH/k < 1 implies 
that / ~ {aH/k)~^aHd{a) > 55 at o ~ 10~^ because 
aHd{a) ~ 55 at last scattering. In summary, for observa- 
tions at large redshifts, the requirement aH/k < 1 provides 
the upper bound on the scales for which our ansatz can be 
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used wheres I > 30 does the same as low redshifts. This up- 
per bound can be relaxed depending on the range of redshift 
in which the observations are made. 

Now, for the lower bound on Ho/k we get Ho/k > 
3 X 10~^. This corresponds to k^i ~ 10~^/iMpc~^ which 
is at the boundary between linear and nonlinear evolu- 
tion of Sm today. At these scales the linear and nonlin- 
ear matter power spectrum differ by a few percent today 
(and less in the past). Since the scalar metric fluctuations 
©[^(x, a), 'I'(x, a)] ~ 10~^ on these scales, as indicated by 
measurements of the cosmic microwave background (CMB), 
we can linearize the field equations in <&, ^ and 5m at these 
scales. Another reason for this lower bound is that on scales 
larger than these we do not expect a significant bias be- 
tween the baryonic and nonbaryonic matter. We can re- 
lax the lower bound if the observations are restricted to 
smaller scale factors since the scale factor dependence of the 
boundary between linear and nonlinear evolution is given by 
fcni(a) ~ W~^a~^^^ hMpc~^. For example if we restrict our 
selves to 10"^ < a < 10"\ then Ho/k > fO~*. 

Figure 1 shows the typical order of magnitude of aH/k 
for the range 3 x 10"^ < Ho/k < 3 x 10"^ (filled yellow 
region). Finally, the range of scale factors we have in mind 
for our framework is lO"'^ ^ ^ 1- Gravitational dynam- 
ics at late times (large a) is particularly interesting due to 
cosmic acceleration. The next generation of observations in- 
cluding lensing, BAO, cluster counts, galaxy power spectra 
etc. will be made within this range. Although we concentrate 
on late times, with some care, our framework can be used 
with CMB observations. For example, after including radia- 
tion and baryons, using our framework we can calculate the 
anisotropics in the CMB if we know the initial conditions 
for each mode after it enters the horizon. Once the modes 
are sufficiently subhorizon, their subsequent evolution can 
be used to constrain the coefficient functions. Note, that for 
the mode corresponding to the first acoustic peak {I ~ 220), 
aH/k ~ 0.3 at last scattering. This comoving scale (as well 
as a range of smaller scales) is within the yellow shaded 
region in Figure 1. 

Before we end this section we provide a concrete exam- 
ple of what the coefficient functions look like in a simple 
case, the Einstein-de Sitter universe; 



(4) 



where ~ fO ^. We turn to the calculation of the coefficient 
functions in the next section. 



3 APPLICATION OF THE FRAMEWORK 
WITH EXAMPLES 

In this section we calculate the coefficient functions for GR 
with a cosmological constant and nonrelativistic matter, GR 
with quintessence, scalar-tensor theories, f{R) theories and 
DGP gravity. In general, the nonrelativistic matter consists 
of baryons, massive neutrinos and nonbaryonic dark matter 



Po 


= 70 = 1, 




= 71 = 0, 


132 


= 72 = -3, 


So 


= a/ai, 


5i 


= 0, 


52 


= 3(a/ai)(l 



with (possibly) nongravitational interactions between them 
and other fields. For simplicity we will ignore massive neutri- 
nos and baryons in this section. Local tests of gravity pro- 
vide strong constraints on baryons and photons and their 
interactions. They do not yet provide similar constraints on 
the interactions of nonbaryonic matter. Hence, nonbaryonic 
matter need not be covariantly conserved. However in the 
examples considered, we treat dark matter as a perfect fluid 
that is covariantly conserved for simplicity. This allows us 
to use the conservation equation ([3|, which is sometimes 
easier to use than a gravitational field equation that would 
otherwise take its place. 

The basic strategy is to substitute our ansatz into the 
field equations and conservation equations and solve for the 
coefficient functions. We begin by substituting our ansatz 
^ into the conservation equation for nonrelativistic perfect 
fluid dark matter (|3|, collecting terms with like powers of 
(aH/k) and setting their coefficient terms equal to zero to 
obtain 

[a^dl + (2 + q)ada] 5o - ^^^/3o<5o = 0, 



[a'^dl + (2 + q)ada] [{aH)5 
-(aH)IJido, 



i/2 



-fio\{aH) 52\ 



^dl + (2 + q)ada] [{aHfS2] 



{aHf[2f3^5i + l325o - ^{a'dt + qada - g)(7o<5o)], 



H^ 



(5) 

where q{a) and H{a) are assumed to be known from the 
background evolution. The above equations are second or- 
der differential equations for 5o, 5i and 52. The equation for 
5o can be solved once /3o is known. G/3o is the effective grav- 
itational "constant". If /3o ~ 1, the equation for 5o is the 
usual equation for the fractional matter overdensity on lin- 
ear and small scales in GR with nonrelativistic matter as 
the only clustering component. 

We digress a bit to note that for 5„ = {aH)"5n, the 
differential operator acting on 5n is [a^da -|- (2 -|- q)ada — 
4-KGpmPo/H'^]. This feature continues if we were to go to 
higher order terms as well, hence it might be useful to find a 
Green's function for this operator. In general, to solve for 5i, 
we need to know Po,l3i and So with two initial conditions. 
Similarly, to solve for 52 we need to know /3o, 70, 5o, 5\ and 
(32 along with two initial conditions. To progress further we 
turn to specific theories of gravitation. Our aim is to show 
how to apply the formalism rather than discuss in detail the 
various models considered. We leave out the detailed steps, 
which are straightforward but tedious. 



3.1 General relativity with cold dark matter and 
the cosmological constant 

We start with the usual Einstein Hilbert action: 

5 = ^ j d\^g [R - 2A] + I d\^£m, (6) 

with Crn, the lagrangian density for perfect fluid cold dark 
matter The corresponding field equations are 



G^+ASi:; = SnGTt:, 



(7) 
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where GfJ = Ri ~ 5'iR/2 and T^' is the energy-momentum 
tensor for a pressureless perfect fluid. As usual, we separate 
the field equations into the background and perturbed parts 
(first order in $ and 5m)- Upon substitution of our ansatz 
into the perturbed field equations we get the following ex- 
pressions/equations for the coefficient functions. 



[a^dl + (2 + q)ada] So - ^^g^5o = 0, 



[a^dl + (2 + q)ada] [{aH)Si] - 
[a^dl + (2 + q)ada] [{aHfSi] - 



[(aH)5^] = 0, 



= jp^ — [aH) [a da + {q+l)ada-q\oo, 

Po = 70 = 1, 
Pi =11= 0, 



P2 



72 



.adaSo 



where we used the 00 and i j Einstein equations along 
with the coefficient form of the conservation equations ([5|. 
We need to provide 6 constants of integration for the three 
second order differential equations. We take these to be 



So{ai) = 1, aidaSo{ai) = 1, 
(5i(a,) = 0, aida6i{ai) = 0, 
<52(t[i) = 0, aida52{ai) = —3. 



(9) 



This ensures that (5m(k, a^) = (5mi(k), thus defining 5mi(k) 
in our ansatz Q. The derivatives are chosen to agree with 
the case of pure matter domination at early times (a^ ~ 
10~^), where the explicit solution takes the form Sq = a/a^, 
5i = and S2 = 2>{a/ai)(l — ajat) after rejecting the de- 
caying modes. For any model under consideration, we can 
choose fix initial condition by rejecting the decaying mode. 
For simplicity, we shall use the above initial conditions for 
the scalar-tensor as well the braneworld models for which we 
plot the coefficient functions. For these models the parame- 
ters have been chosen so that at ~ 10~^, the conservation 
equations approach those of an Einstein-deSitter universe in 
GR. 

The dashed lines in Figures 2 and 3 show these di- 
mensionless coefficient functions for the spatially flat-ACDM 
with Q.ra = SivGpmo/SHQ — 0.3. Since Po = Jo = I, there 
are no corrections to the Newtonian gravitational constant 
as far as growth of perturbations is concerned on small 
scales. Since single gradients do not appear in the Einstein 
equations involving Sm , "1> and 9 (after eliminating the ve- 
locity through the conservation equation). Pi = "fi = 0. 
The 00 Einstein equation imposes Si = 0. The fact that 
/32 = 72 7^ reflects corrections because of GR to the rela- 
tionship between matter and metric perturbations, whereas 
/32 = 72 7^ —3 reflects the effect of the cosmological constant. 
So characterizes the growth of structure on small scales. It 
deviates from Sq = a /at because of A. 82 reflects the cor- 
rections to the growth function as we move to larger scales. 
Note that P2 and 72 and 82 are multiplied by (aH/k)^, whose 
magnitude is shown in Figure 1. The terms P2{aH/k)^, 
12(0,11 / k)'^ and S2{aH/k)^ are much smaller than Po, 70 and 
So, making it difficult to observe their effects unless we in- 
vestigate large scales. 



3.2 Scalar-tensor theory with cold dark matter 
(matter representation) 

Scalar-tensor theories are popular alternatives to GR. In the 
matter representation (also called the Jordan frame), the 
action contains two free functions /((^) and V{ip) 



S = 



(10) 



Note that we have decided to make dimensionless since 
we wish to treat the perturbation in this field Sip on the 
same footing as the metric perturbations "I> and ^. Also 
C^p = —{8^)^/2 — V{tp) and Cm does not contain The 
field equations for this theory are 



8^G 1 



d^ipdvp — Sii 



f 2/ 
The field equation for ip is 

n<p-v^ + UR = 0, 



-d^ipdap + V 



(11) 



(12) 



where = d^f and 1/^ = di^V . 

These field equations at the background level can be 



found in the literature (for example see Boisseau et al. 2000 1 



Using our ansatz in the perturbed gravitational field equa- 
tions and the field equations for ip at first order in $, 5m 
and Sip>, collecting terms with like powers of (aH/k), and set- 
ting the expression in front of each power of {aH/k) equal 
to zero, we get the following expressions/equations for the 
coefficient functions; 



1 /1 + 4/a^ 



Po = 



Pi=li= 0, 



/ V 1 + 3/a2 
1 fl + 2fa^ 



+ 0[c 



P2 



3 adaSo , 1 . „ \2 

j^ + ^(aa.^) 



+ 



-2,{ada^) 



adaSo 



- {adapY ^ + Ka-da^p) + ^ 



72 



+ 0[a% 
3 adaSo 

'^l~sr 



, . adaSo 1 . a \2^v 
00 2 a 



+ 



{ada<p) 



2H^ 



a 

7 

(13) 



where a = f^p/ f is the coupling function and all the func- 
tions depend on the scale factor a . We have calculated the 
full expressions for P2 and 72, which are rather long. The 
first two terms are listed as a power series in the coupling 
function a <C 1 with We used the i / j 

equation, aSip = — to eliminate Sip from the held equa- 
tions. The 00 equation and the held equation for Sip yield 
P„ and 7„ with (n = 0,1,2). The equations for So, Si and 
S2 are given by equations ([5| with P„ and !„ (n — 0, 1,2) 
given above. Again using the initial conditions ([9|, we can 
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0.67 0.25 



0.01 0.2 



Figure 2. The dimensionless coefBcient functions characteriz- 
ing the relationship between the metric perturbations and matter 
distribution are shown above for i<'ACDM(dashed Hues) and the 
scalar-tensor theory (STT) (solid lines). The STT model is cho- 
sen so that its expansion history is consistent with observations. 
In the case of ACDM /Jg = 70 = 1, /3i = 71 = and f^2 = 72- 
At early time (matter domination) /32 = 72 = ~3 with the cos- 
mological constant causing a departure from this value at late 
times. The variation of fSo with the scale factor in the STT can 
be interpreted as a variation of Newton's constant "G/3o" as far 
as growth of perturbations is concerned. Also note that for STT, 
00 7^ 70 S'lid 02 7^ 72- For STT, the difference in the coefficient 
functions is due to $ — 4* = —a{ifi)Sip ^ 0. Note that /3i = 71 = 
in STT as well as ACDM. We remind the reader that in the ansatz 
^ the coefficients (^2 and 72 are multiplied by (aH/k)^, whose 
magnitude is shown in Figure 1, making them accessible at large 
scales only. 



120 




a 



Figure 3. The dimensionless coefficient functions characterizing 
growth of structure are show above for ACDM(dashed lines) and 
the scalar-tensor theory (STT) (solid lines). The STT model is 
chosen so that its expansion history is consistent with observa- 
tions. So is the usual growth function on small scales, whereas 
62 characterizes the departures as we move to larger scales. For 
ACDM and STT, 5i = 0. We note that S2 is the coefficient of 
(aH/k)^, which is small withing the scales of interest (see Figure 
1). The initial conditions for Sq and 62 are chosen at ai ~ 10~^ 
and are consistent with growth of structure in a matter dominated 
era. 



3.3 General relativity with cold dark matter and 
quintessence 

GR with quintessence is a special case of the scalar-tensor 
theories discussed above, with /(<^) = 1- The action and 
corresponding field equations are 



167rG 



(14) 



solve for all the coefficient functions once fif) and V{ip) 
have been provided. Note that the difference <1> — \l/ depends 
on l3n — "fn {n — 0, 2). This is usually small for q ^ 1 since 
/9o — 70 ~ and /32 — 72 ~ a- 

We plot the coefficient functions in Figures 2 and 3. 
We have chosen f{ip) = 1-1- ci(p'^ and V{ip) = 2A(1 -I- C2(p'^) 
with ci = C2 =0.1. The initial conditions and ci,C2 were 
chosen to ensure that the expansion history remains consis- 
tent with observations (consistent with ACDM to within a 
few percent). The difference between /3„ and 7„ (n = 0,2) 
is due to nonminimal coupling {a 7^ 0). We stress that we 
have not included baryons in this illustrative calculation. In- 
cluding baryons would lead to very strong constraints on the 
function f{ip) today from solar system tests ( [Schimd et al.| 
2005j). For an example of a STT that includes dark matter 



and baryons with different couplings to gravity see (Bean 
[2001] ). 



0^ = 8nGT^ + 



□V5 - = 

The coefficient functions are given by 



(15) 



(16) 



l3o 
I3i 



70 = 1, 

71 =0, 

72 



adaSo 1. T s2 

do 4 



(17) 



where {ada(p)^/4: = 1 — q — AivGpm/H'^- The i / j Einstein 
equation yields /3„ = 7n (n = 0, 1, 2). We used the Oi equa- 
tion to eliminate Sip from the field equations. As before So, Si 
and ^2 are provided by equation 

We pause to comment on a difference between min- 
imally and nonminimally coupled scalar-tensor theories. 
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Consider the field equation ( 12 1 for Sip 

\2 



[a^dl + (3 + (j)a9„l Sip + 



k 
oH 



+ 



6(1 + q)U 



Sip 



{ada<p ~ 6fa)ada^ + S{adaP - 2(4 + q)fa)ada'f 



2 (^6/a(l + g) + J $ + 2/a (^^J ($ - 2*) 

(18) 

In the minimally coupled case we set /{f) = 1, a{ip) = to 
get 

\2 



[a'^dl + (3 + q)ada] Sip + 
^{ada(p)ada^ + 3{ada(p)ada'^ 



k 



+ 



VP 



Sp 



(19) 



From the above equations we can see that in the nonmin- 
imally coupled case, for k/aH ^ 1 we have 5p oc a{p)^ 
whereas in the minimally coupled case 5p cc ${aH/k)^ . 
Along with $ oc Sm{aH/k)^ , at large k the additional field 
5p follows the same aH/k expansion as the potentials with 
Smi(k) multiplying the expansion. This is one of the argu- 
ments we had used in Section 2 to justify the form of our 
ansatz. We have assumed 5' = 0[$] in this argument. 



3.4 f{R) gravity with cold dark matter 

In recent years modifications of the Einstein-Hilbert action 
in the form of a function of the Ricci scalar has become a 



popular alternative to quintessence (see for example Carroll 



et al.||2004[ [Nojiri fc Qdintsov||2007| ). The action and field 
equations are 



S 



16tvG 
1 + fR)G[ 



d\^[R + fiR)] + 



(20) 



In the above expressions /_r = 9j;/(i?). The coefficient func- 
tions are 



l3o 



3(1 + /fl) 
2 



3(l + /fl)' 
/3i = 71 = 0, 



(1 + Jr) 



2 a^dlSo 

3 So 



72 



-lmB{j + q-2) + 2-q}'^ 
6 Oq 

{10 - 4g + + 2j(g - 4) - s} 
+72B\j + g - 2)' - 4(j + g - 2)adaB + 2q 
1 



1 

9B 



(l + /«) 



2a^diSo 
3 So 



T-dgSo 
So 



+ ^{{6B(j+q-2)-7 

-4B{4-g + g2+j(2g-5)-s} 

1 



+4:{j + q- 2)adaB - 2q + 



9B 



(21) 

where j — dg/dlna— (1— 2g)g and s — dj/d\na—(2—3q)j are 
the scale factor dependent functions, jerk and snap respec- 
tively, and B = H^fRR/{l + To obtain So,5\ and ^2 
we use equation ([5|. Again, as in the case of GR and scalar- 
tensor example, the coefficients of the aH/k term in the 
ansatz vanish. Note that we have assumed (aH/k)^ B^^ ^ 1 
in deriving the above expressions, hence it is not appropriate 
to take the limit i? — > after deriving the coefficient func- 
tions. Under this assumption, to lowest order in aH/k, we 
get <E> = 2*, unhke GR with ACDM where $ = *. If we take 
the opposite limit, (aH/k)^ B~^ ^ 1 the coefficient func- 
tions are quite different. In particular, /3o = 7o = (1 + Jr)^^ 
and we reach the GR limit as we let fR—>0. As long as we 
ensure, a priori, that this transition scale (H/kc) ~ B^^^ 
(see for example Pogosian & Silvestri]'2008 ) is outside the 
length range of physical scales of interest, wc can use our 
ansatz. More details on the dynamics of f{R) theories in 
the context of structure formation, solar system tests, etc. 



d in ( Faulkner et al. 


2007 Songet al.|2007 Bean 


|Chiba et al.||2007| 


Pogosian & Silvestri 2008|). 



Finally we note, that our purpose in discussing f{R) mod- 
els was to illustrate an application of our framework. These 
models suffer from a number of problems including fine tun- 
ing to match the solar system constraints as well as a rather 
serious instability, where the curvature blows up at finite 
matter densities ( |Frolov|2008[ ). 

3.5 Brane world models: DGP Gravity 

As a final example, we provide the expressions and equations 
governing the coefficient functions for DGP gravity. In this 
model, matter is restricted to a four dimensional brane in a 
five dimensional bulk. In addition to the the Einstein-Hilbert 
action in the bulk, there is an induced four dimensional term 



( |Dvali et al. 2000 1. More explicitly, the full 5D action is given 
by 



32nGr^ 
1 



+ 



167rG 



J d^x^-g^^s)R(&) 
^xy/—gR 



(22) 



In the above action rc ~ G'^^^ /2G where G'^' is 5D gravita- 
tional constant. The field equations are given by the Einstein 
equations in the bulk (A, B = 0, 1, 2, 3, 4): 



Gi 



0, 



(23) 



and the Israel junction conditions on the brane (/i, v 
0,1,2,3) 



-QS'i 
3 



(24) 



where Ki^ is the extrinsic curvature on the brane. In Gaus- 
sian normal co-ordinates, the extrinsic curvature is given by 
the derivative normal to the brane 



(25) 



^ Our B = H'^ fnn/ (1 + fu) differs from tfio definition of B in 
l |Song et al'!]|2007| by a factor of (q - l)/6(j + q - 2). Also it is 



unrelated to the B used in 
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Figure 4. The dimensionless coefficient functions characteriz- 
ing the relationship between the metric perturbations and matter 
distribution are show above for ACDM(dashed lines) and DGP 
braneworld model (solid lines) . The variation of /3o with the scale 
factor in DGP can be interpreted as a variation of Newton's con- 
stant "G/3o" as far as growth of perturbations is concerned. Also 
note that for DGP, /3o 7^ 7o • In contrast to all the other examples 
considered, the coefficients of aH/k: 71 ^ 0. This is due to the 
junction conditions on the brane. The linear aH/k term provides 
an intriguing signature of braneworld models. 



On the RHS of the junction conditions 



(26) 



where GJJ and are the 4D Einstein and stress-energy 
tensors respectively. For this model our coefficient functions 
are given by (which can be easily determined from the results 
in |Deffayet|2002l ) 



l3o 



70 



Pi = 



71 



A-2Hrci2 + q) 
3-2Hrc{2 + q)' 
2 - 2Hr,{2 + q) 
3-2Hrc(2 + q)' 
12(-1 + Hr,{l- 



1)? 



(27) 



Hrci3~2Hrc{2 + q))2' 
6{1 - 2Hrc){l - Hr^il + q)) 



Hr^{3-2Hr,{2 + q)Y ' 

We can solve for So and 5\ using (|5|. Note that even though 
5i{a.i) = aidaSi{ai) — 0, 5i{a) 7^ because /3i / 0. We note 
an important difference between the the DGP braneworld 
model and the examples considered so far in this paper. 
Unlike the previous examples, the coefficients of the odd 
power of aH/k are non-zero (/3i,7i,5i 7^ 0). As explained 
below, the odd power of aH/k arises due to the junction 
conditions that must be satisfied by metric perturbations 
at the location of our four dimensional brane in the higher 
dimensional bulk. 

A general way to understand the odd power in our 
{aH/k) expansion is as follows. The Israel junction condition 
relates the first (y) derivative of a metric perturbation nor- 




Figure 5. The dimensionless coefficient functions characterizing 
growth of structure are show above for ACDM(dashed lines) and 
DGP braneworld model (solid lines). Sg is the usual growth func- 
tion on small scales, whereas 5i characterizes the departures as 
we move to larger scales. In contrast to ACDM, for the DGP case 
(5i 7^ 0. This could provide a distinct signature of braneworld 
models. 



mal to the brane at its surface to the 4D Einstein tensor and 



stress-energy tensor in the brane (see equations (24 1,(25 1 
and (26l). The 5D vacuum Einstein equations in the bulk 



provide homogeneous second-order linear differential equa- 
tions for the metric perturbations. Just outside the brane, 
the operators in these equations will be dominated (in our 
large k limit) by the 3-dimensional spatial derivatives (con- 
tinuous across the brane surface), giving terms proportional 
to fc^. These must be balanced by terms proportional to dy. 
Thus we see that dy must be proportional to k. 

We note that the above arguments are rather general. 



Although the RHS of equation ( 24 \ will be different in dif- 
ferent braneworld models, due to the LHS a linear k~^ term 
will be present in most braneworld models. However, the co- 
efficients might not be of 0[1] as in the case of DGP. The 
existence of this non-zero linear aH/k term provides an ex- 
citing new signature for the DGP (or other ) braneworld 
models. Since on subhorizon scales aH/k ^ 1, it is signifi- 
cantly easier to constrain the linear term compared to the 
quadratic one. 

We plot the coefficient functions in Figures 4 and 5. For 
DGP, we have assumed flm = 0.25 and Qt ~ for simplicity 



although this model is in tension with current data (Song 
et al. 2007[ |. The dotted lines represent the coefficients in 
ACDM. We have not calculated the coefficients /32,72 and 
52 of the [aH/kf terms for the DGP case because they are 
expected to be subdominant compared to the linear aH/k 
terms. 



^ Another way of seeing the odd power of aH/k in DGP case is 
through the form of a propagator which involves ( Dvali et al. 



2000 1 
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We have ignored two important features in the DGP 
model, the strong couphng problem and the ghost prob- 



lem in the self accelerating branch (see for example Deffayet 
et al.||2002[|Luty et al.||2003[ ). The transition to the strong 
coupling regime happens at the Vainshtein radius r* . For a 
localized matter distributions with Schwarzchild radius r 



^2x1/3 



Using 



1/Ho, for the largest localized 



distribution in our Universe, galaxy clusters {M ~ 10 Mq), 
we get < lOMpc. This scale is well below the smallest 
scales where we intend to apply our framework. We are un- 
aware of a calculation for r^, when considering distribution 
of matter on cosmological scales (which is not localized) . To 
avoid the ghost problem one can choose the normal branch 
rather than the self-accelerating branch of the DGP model 
(see for example Charmousis et al.|2006 1 . A calculation sim- 
ilar to the one done in this section can be repeated for the 
normal branch, however in that case we do not have an ac- 
celerating universe. 

In this section we have calculated the coefficient func- 
tions for a few examples. Our aim was to give a flavor of the 
calculations rather than be exhaustive in the investigation of 
the models considered. It would be interesting to investigate 
these models in more detail in the context of these coefficient 
functions to see if there is come generic behavior across a 
large class of models. Based on the examples considered it 
might be tempting to conclude that /?„ — 7n 7^ indicates 
physics beyond general relativity. However this is not so. For 
example a hypothetical dark energy component could also 
yield significant anisotropic stress. In the early Universe, a 
more standard source of anisotropic stress was provided by 
neutrinos. Nevertheless this difference could serve as an indi- 
cator of new physics in the matter or gravity sector. Another 
intriguing result was the presence of a term linear in aH/k 
in the braneworld scenario, which could provide a unique 
signature of braneworld models. We have left out many pos- 



sibilities including Bekenstein's TeVeS (Bekenstein 



models with non-canonical kinetic terms ( Chiba et al 



20041, 



20001, 



models of imperfect fluid dark energy with anisotropic stress 



(Koivisto & Mota 20061, and many others (see Copeland 
et al.|2006 for a review). We now turn our attention to ob- 
servables and their relationship to the coefficient functions. 



4 OBSERVATIONAL IMPLICATIONS 

We have outlined a procedure that allows many alternative, 
dynamical theories to GR with FACDM cosmology to be ex- 
plored within a common framework. Our approach has been 
devised with future observations in mind as its usefulness is 
limited to the observations that we expect will be the most 
prescriptive. We assume that the expansion history of the 
Universe is well constrained through the distance redshift re- 
lation obtained from the apparent magnitude of Type la su- 
pernovae, the baryon acoustic oscillation scale and the ratio 
of baryons to total matter in galaxy clusters. The large k ex- 
pansion connects the inhomogeneous nonrelativstic matter 
distribution to the perturbed metric in a universe of known 
(unperturbed) kinematical behavior, i.e. with a given rela- 
tion H{a) [or, equivalently, a{t)]. Our approach also pre- 
sumes that the theories under consideration provide an un- 
derstanding of how the distribution of observable entities 
such as galaxies relate to that of total mass. This allows us 



to focus on the manner in which structure can be observed 
to grow in the linear regime well within the horizon, which 
avoids the limitations imposed by cosmic variance consider- 
ations and the complications associated with gas dynamics. 
We further suppose that gravitational motion of baryonic 
matter and photons follows timelike and null geodesies re- 
spectively in this spacetime. 

From an observational standpoint, our focus is on co- 
moving length scales from ~ 40 Mpc to ~ 400 Mpc or equiv- 
alently 300 > i > 30 at 2 ~ 1, where we expect the effects 
to stand out the best. There are three types of observations 
that are likely to be relevant. Firstly, there are direct mea- 
surements of the two point correlation function and its evolu- 
tion. Counting galaxies (or clusters) in three dimensions will 
lead to measurements of the evolution of the density function 
5m using future survey instruments such as LSST (Tyson 



|2002{ |Zhan et al.||2006| ) limited solely by cosmic variance as 
the photometric redshift accuracy and biasing errrors will 
be ignorable on these scales. We can construct the ratio of 
the matter power spectrum P^^ {k, a) at different redshifts 
to obtain constraints on 5o{a), 5i{a) and S2{a). As discussed 
in the Introduction, by taking ratios, we can eliminate the 
need for knowing the initial conidtions (5mi(k). 



Ps^{k,a2) 



5o{ai) 



where {/(fc, a)}^J = f{k, 02) - /(fc, ai). 

The second type of observation that will be carried out 
involves departures from the Hubble flow. These are domi- 
nated by the potential function $. Under our assumptions, 
galaxies will follow timelike geodesies and satisfy the linear 
conservation equations relating their peculiar velocities to 
<!>. 

Finally there are weak lensing observations which de- 
pend upon the sum, $ -I- ^P, presuming photons follow null 
geodesies. These then allow us to track the evolution of A 
combination of these measurements would not only allow us 
to understand the scale dependent evolution of $, 5' and S,n 
but also allow us to probe the relationship between them. 
For example, using our ansatz, one can obtain constraints 
on the coefficient functions by comparing the correlation 
functions for the potentials, P$-|-*(A;,a) (provided by lens- 
ing tomography) and the nonrelativistic matter overdensity 
Ps^ {k, a) (provided by growth of structure measurements) 
using 



fe*P$+4, oc Ps^ {/3o +70)^ 



1 + 2 



Pi +-71 



aH 

IT 



, A) + 70 , 

Comparing the matter and potential power spectrum allows 
us to constrain the coefficient function without worrying 
about the initial conditions, though one would still have to 
obtain this ratio at different redshifts to constrain the time 
evolution of the coefficient functions. 

In this exploratory paper, we have discussed only a 
handful of observations that can allow is to constrain the co- 
efficient functions. In addition to the observations mentioned 
above, we list a few other observations that we think might 
be relevant for our framework. The matter and potential 
fluctuations at the last scattering surface can be compared 
to their counterparts at late times, as long as we restrict our- 
selves to linear subhorizon scales. The same is true for BAO 
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measurements (see discussion of range of scales at the end of 
Section 2). Recently, a 3cr detection of lensing of the CMB at 
large I, has been reported by the ACBAR group (Reichardt 
et al.|[2008 K This measurement probes the distribution and 
evolution of potentials after last scattering, and can also 



be used for constraining the coefficient functions ( Calabrese 
et al.|2008[ ). With the Planck mission (jPlanck-CoUaborationl 
2006[ ), such constraints are expected to improve significantly. 
Another exciting probe of the three dimensional matter dis- 
tribution may be provided by the 21 cm observations (see 
for example |Mao et al. [2008] and references therein). 

We have limited ourselves to the linear regime. On small 
scales, the nonlinear matter power spectrum and its evolu- 
tion can play a role in the observations discussed above. The 
linear to nonlinear mapping discussed in ( Smith et al.|2003 1 
can be used for this purpose. Ifowever, without understand- 
ing the theories under consideration in the nonlinear regime, 
this is not fully robust. 

Recall that {/9n,7n,<5n} with (n = 0, 1,2) are functions 
of the scale factor, a. If the observations are to be done in 
a limited range of redshifts then Taylor expanding the co- 
efficient functions around the central value of the redshift 
might be a simple and model independent way of character- 
izing these coefficient functions in terms of a few parameters. 
From a theoretical perspective, the coefficient functions will 
depend on relevant parameters in the theory or model under 
consideration. A detailed investigation of the parameteriza- 
tion of the coefficient functions and the possible constraints 
that can be obtained from current and future observations is 
beyond the scope of this paper. For a more detailed discus- 
sion of the observations for distinguishing different models 
of modified gravity and dark energy we refer the reader to 
some of the references cited at the end of Section 1 in this 
paper. 



5 DISCUSSION 

We have outlined a procedure that can be used to test the 
application of general relativity (more specifically FACDM) 
on cosmological scales in the context where it is most likely 
to fail and in the regime where observations should be most 
sensitive to measuring a departure from the general rela- 
tivistic prediction. The scales are large enough to avoid the 
complications from nonlinearities and gas physics, yet small 
enough to avoid strong limitations to the interpretation of 
observations posed by cosmic variance. 

Our procedure assumes that (i) The geometry and kine- 
matics of the Universe is understood (ii) baryons and pho- 
tons behave as ideal test particles following geodesies of 
the cosmological metric. Given these assumptions, at late 
times, it is the relationship between the cosmological metric 
and the nonrelativistic matter distribution (along with their 
respective evolution) that provides a test for alternatives 
GR with a cosmological constant and cold dark matter. To 
probe the dynamics of gravity (or any additional fields) we 
provided an ansatz, equation(|2|, which gave a relationship 
between the cosmological metric and nonrelativistic matter 
perturbations in the linear, subhorizon regime. This form of 
the ansatz is consistent with a large class of theories with 
the differences between different theories evident in the co- 
efficient functions {/3„(a), 7„(a), (5„(a)} with n = 0,1,2. ft is 



hoped that three scalar functions, the nonrelativistic matter 
overdensity Sm and the metric potentials $ and ^ can be 
measured over the next decade, providing constraints on the 
coefficient functions. Constraining these coefficient functions 
provides observers with concrete targets for testing gravity 
in a scale dependent manner. 

Our goal was to provide a perturbative framework, sim- 
ilar in spirit to the PPN formalism for testing gravity on 
solar system scales. However unlike the PPN case, we were 
left with coefficient functions that depend on the scale factor 
rather than constant coefficients. Although we have not done 
so in this paper, if the observations are limited to a small 
range of scale factors, it is possible to characterize these 
coefficient functions using a few parameters by expanding 
around a given scale factor at which the observations are 
centered. 

With our choice of scales, we have restricted ourselves 
to linear, subhorizon evolution. We leave the connection be- 
tween superhorizon and subhorizon evolution as well as con- 
sideration of nonlinearities for the future. Although, we have 
restricted ourselves to scalar perturbations, the framework 
could be extended to include vector and tensor perturba- 
tions. 
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